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ABSTRACT 

Having already reported on the first photometric study of the intermediate-age open 
cluster NGC 1883 (Carraro et al. 2003), we present in this paper the first spectroscopic 
multi-epoch investigation of a sample of evolved stars in the same cluster. The aim is 
to derive the cluster membership, velocity and metallicity, and discuss recent claims 
in the literature (Tadross 2005) that NGC 1883 is as metal poor as globular clusters 
in the Halo. Besides, being one of the few outer Galactic disk intermediate-age open 
clusters known so far, it is an ideal target to improve our knowledge of the Galactic 
disk radial abundance gradient, that is a basic ingredient for any chemical evolution 
model of the Milky Way. The new data we obtained allow us to put NGC 1883's 
basic parameters more reliable. We find that the cluster has a mean metallicity of 
[Fe/H] = -0.20±0.22, from which we infer an age (eSOly" Myr) close to the Hyades 

one and a Galactocentric distance of 12.31q'2 ^P<^- The metal abundance, kinematics, 
and position make NGC 1883 a genuine outer disk intermediate-age open cluster. We 
confirm that in the outer Galactic disk the abundance gradient is shallower than in 
the solar vicinity. 

Key words: Open clusters and associations: general- Open clusters and associations: 
individual: NGC1883- Milky Way: structure and evolution 



1 INTRODUCTION 

NGC 1883 is a northern open cluster, located in the sec- 
ond Galactic quadrant, toward the anticenter direction (a = 
05''25'".9, 5 = -H46''29', I = 163°.08, b = -h06°.16, J2000.0). 
This cluster was studied by Carraro et al. (2003, hereafter 
C03), who present the first BVI CCD photometric observa- 
tions suggesting an age of about 1 Gyr and a distance of 4.8 
kpc from the Sun. 

This would imply a Galactocentric distance larger than 13 
kpc, and would make NGC 1883 a key object four our under- 
standing of the properties of the outer Galactic disk. How- 
ever, both age and distance are uncertain because of the 
unknown metallicity of the cluster. 

At the position of NGC 1883 the mean abundance (Friel 
1995) should be significantly lower than solar. Unfortu- 
nately, the exact amount of this under-abundance is poorly 
constrained, also in light of the results by Twarog et al. 
(1997), Carraro et al. (2004), and Yong et al. (2005), that 
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suggest that the disk abundance gradient for Galactocentric 
distances larger than '^lO kpc gets much shallower than in 
the solar vicinity. 

Using solar metallicity (Z=0.019) C03 estimated E(B- 
V)=0.23, (m-M)y=14.00 for NGC 1883, and E(B-V)=0.35, 
(m-M)y = 14.50 adopting a lower metal content (Z=0.008), 
which by the way seems to provide a better fit to the star 
distribution in the Colour Magnitude Diagram (CMD). 
A more recent study was performed by Tadross (2005, here- 
after T05) using C03 BVI data and JHK 2MASS data. 
This author obtained a photometric metallicity of about 
[Fe/H]~-l.l (Z~0.0015), much lower than that estimated 
by C03, and closer to globular clusters. The disagreement 
between the two abundance estimations has to be solved in 
order to obtain reliable parameters for this cluster. In fact, 
NGC 1883 is one of the few intermediate-age open clusters 
located in the outer Galactic disk. For this reason, it plays 
an important role in defining the shape of the radial abun- 
dance gradient in the external regions of the Galactic disk, 
which is one of the fundamental costrains of models of chem- 
ical evolution (Cescutti et al. 2007). 

For this purpose, this paper presents the first spectroscopic 
study of the cluster aiming at finding the member stars 
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from radial velocity measurements, and to estimate a re- 
liable value for the mctallicity. Section 2 of this paper illus- 
trates the observation and reduction strategies. In Section 3 
and 4 we discuss the membership and the abundance deter- 
mination. Section 5 is devoted to a comparison between this 
paper and previous ones. Section 6 deals with a new deter- 
mination of cluster reddening, distance and age. Section 7 
compares the observed CMD with a synthetic one. Section 
8 and 9 are devoted to the determination of structural pa- 
rameters for the cluster (radial profile, luminosity function, 
and mass function). Finally Section 10 discusses NGC1833 
in the context of the Galactic abundance gradient. 



2 OBSERVATIONS AND DATA REDUCTION 

Medium-high resolution {R w 20000) spectra of 5 candi- 
date evolved stars in the field of NCG1883 have been ob- 
tained using the REOSC Echelle Spectrograph on board 
of the 1.82 m telescope of Asiago Astronomical Observa- 
tory. The Echelle spectrograph works with a Thomson 1024 
X 1024 CCD and the wavelength coverage for our obser- 
vations is approximately 4600-6400 A. Details on this in- 
strument are given in the Asiago Observatory Home page 
(www.oapd.inaf.it /asiago/) . 

The exposures were 3x15 min for all stars. To improve the 
signal to noise ratio the three exposures were added reaching 
S/N values up to 20. Each star was observed during three 
epochs, December 2004 (A), January 2005 (B), January 2006 
(C), allowing us to establish the binarity of each target. 

Data were reduced with the IRAF package 1^, including 
bias subtraction, flat-field correction, frame combination, ex- 
traction of spectral orders, wavelength calibration, sky sub- 
traction, and spectral normalization. By comparing known 
sky lines positions along the spectra we obtained an error in 
wavelength calibration of less than 0.01 A (<0.5 km/s). 
Details of the observed stars are listed in Table 1, while the 
positions of these targets on CMDs axe plotted in Fig. 1. 
The first column of Tab. 1 gives the ID number of the star 
according to USNO-2.0 catalogue, the second and third the 
coordinates, the fourth, fifth, sixth, and seventh the V mag- 
nitude and B-V,V-I,V-K colours, the eight, ninth, and tenth 
columns the measured radial velocity at different epochs ex- 
pressed in km/s. 

The photometry for our stars was taken from C03 (B,V,I 
magnitudes), and from 2MASS (K magnitude). 

Radial velocity was measured using the IRAF package 
FXCOR, which cross-correlates the observed spectrum with 
a template. As a template wc used a synthetic spectrum 
having mean atmospheric parameters close to the cluster 
RGB stars (Teff=4500 K, log(g)=3.0, [Fe/H]=0.0, vt=1.5 
km/s). Finally the heliocentric correction was applied. The 
typical error for our measurements is less than 1 km/s. 



t IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 




Figure 1. CMD for the stars in the field of NGC 188.3. left: V vs 
B-V diagram from 003. center: V vs V-I diagram from O03. right: 
V vs V-K diagram (K magnitude from 2M ASS) . Open squares are 
member stars (#75374,#76058), while open triangles are back- 
ground objects. See text for more details. 



3 MEMBERSHIP AND CLUSTER MEAN 
RADIAL VELOCITY 

We make use of these radial velocities to establish star mem- 
bership to NGC 1883. Firstly, we notice that each star has 
radial velocities constant among the three epochs within the 
errors. We are led to the conclusion that none of the tar- 
gets arc short-period binaries. The mean radial velocity and 
membership of all the targets is reported in Tab. 2 where 
member stars are indicated by M, while non-member by NM. 
Only #75374 and #76058 stars have compatible radial ve- 
locity and therefore we consider them as genuine members 
(open squares in Fig. 1) of the cluster. The average of the 
velocities of #75374 and #76058 gives us the mean radial 
velocity of the cluster: 

< RVhelio >= -30.8 ± 0.6 km/s 

The members stars will be used in the following Section to 
obtain the metal content of NGC 1883. 



4 ABUNDANCE MEASUREMENTS 

The [Fe/H] content was obtained from summed spectra of 
member stars (#75374 and #76058) in a interactive way. 
The most important parameter for abundance determina- 
tion is Toff , that for our stars can be obtained from dcred- 
dened B-V, V-I, and V-K colours. But the reddenings E(B- 
V), E(V-I), and E(V-K) for the cluster can be obtained 
only by isochrones fitting, with isochrones having the cor- 
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Table 1. Observed stars 



ID (U- 1350-057...) R.A.(J2000.0) 


Decl.(J2000.0) 


V 


B-V 


V-I 


V-K 


RVh(A) 


RVh(B) 


RVh(C) 


79024 


05:26:02.01 


+46:29:26.9 


12.30 


1.85 


1.89 


4.32 


34.4 


33.0 


34.4 


75374 


05:25:51.69 


+46:29:35.4 


13.25 


1.78 


1.79 


4.00 


-29.7 


-32.5 


-32.1 


76058 


05:25:53.66 


+46:26:20.5 


13.30 


1.61 


1.70 


3.70 


-29.8 


-30.3 


-30.7 


81797 


05:26:09.92 


+46:25:55.4 


11.85 


1.87 


1.94 


4.38 


2.7 


3.0 


2.7 


85248 


05:26:19.92 


+46:31:41.4 


12.61 


2.18 


2.60 


5.43 


-17.9 


-18.5 


-18.4 



rect metallicity. For this reason wc firstly assumed a solar 
metallicity for our stars, obtaining a first guess for reddening 
by isochrones fitting. Then, using the Tefi value calculated 
for this reddening, we measured the metallicity for our stars 
(see below). Using this metallicity, we obtained a new red- 
dening and a new abundance. This procedure was iterated 
until convergence was achieved. 

Teff was derived from the colour- [Fe/H]-temperature re- 
lation from Alonso et al. (1999) using dereddened colours. 
The gravity was derived from the value given by the 
isochrones at the given position of the stars in the CMD, 
while the microturbulent velocity come from the relation 
(Houdashelt et al. 2000): 

vt = 2.22 - 0.322 x log(g) 

The adopted atmospheric parameters for our stars axe re- 
ported in Table 2, and the typical random errors are 100 K 
for Toff, 0.2 for log(g), and 0.2 km/s for vt. 
We derived metallicity from equivalent widths of selected 
spectral lines. The equivalent width was derived by gaus- 
sian fitting of these spectral features. Repeated measure- 
ments show a typical error of about 10 mA. Because of the 
low S/N and the resolution of our spectra, we could measure 
only the strongest Fe lines (~10 for each star) in the region 
A > 5500 to avoid crowding. For this reason we give only Fe 
content of our stars (sec Tab. 2). The LTE abundance pro- 
gram MOOG (freely distributed by C. Sneden, University of 
Texas, Austin; see http://verdi.as.utexas.edu/) was used to 
determine the metal abundances. Model atmospheres were 
interpolated from the grid of Kurucz (1992) models by us- 
ing the values of Teff and log(g) of Tab. 2. The error on 
measurements is high (0.2-0.4 dex), but a weighted mean of 
the obtained metallicities gives a good estimation of the iron 
content: 



[Fe/H] = -0.20 ± 0.22 

This value is important for the discussion reported in the fol- 
lowing Sections and for the comparison with previous works. 
The errors for iron content we report here and in Tab. 2 
represent only internal random errors, mainly due to the un- 
certainty in the equivalent width measurements. Systematic 
errors are more difficult to estimate, and are mainly due to 
uncertainties in colour- Tcff relations and in reddening. The 
typical systematic error on Tefi scale for this temperature 
range is 100 K giving a systematic uncertainty on metallicity 
of 0.02 dex. 



5 COMPARISON WITH PREVIOUS WORK 

Firstly, we compare our results with COS. In that paper the 
authors estimate a range in metallicity between Z~0.019 
([Fe/H]=0.0) and Z~0.008 ([Fe/H]=-0.41) by isochrones fit- 
ting of the TO. Our measurements lie exactly in the mid- 
dle confirming those results. On the other hand, our value 
largely disagrees with T05. This author gives [Fe/H]~-l.l, 
based on the slope of the RGB stars (Vallenari et al. 2000 
and references therein). The method can give a reasonable 
estimate of the metallicity of a cluster (Tiede et al. 1997), 
but [Fe/H]~-l.l is definitely ruled out by our measurements. 
T05 results can be explained by a poor treatment of stellar 
contamination. As a matter of fact, background stars always 
contaminate the field of a cluster, and an identification of 
the members is needed before applying RGB-slope method. 
This is the case of NGC 1883, where a fraction of stars that 
apparently populate the RGB are background objects, as 
Fig. 1 unambiguously shows. This contamination changes 
the real slope of the RGB and the [Fe/H] estimate by T05 
is clearly affected by this problem. As shown in Fig. 1 the 
RGB slope in T05 is mainly based on #79024, #81797, and 
#85248 .star.s (open triangles) that, according to our radial 
velocities, are not members. This leads us to the conclusion 
that T05 metallicity is untenable. 



6 CLUSTER FUNDAMENTAL PARAMETERS 

Having an estimate of the metal content ([Fe/H]=-0.20), we 
now provide a more reliable estimate of the cluster parame- 
ters. We make use of the comparison between the distribu- 
tion of the stars in the CMDs presented in C03 and a set 
of theoretical isochrones from the Padova group (Girardi 
et al. 2000). In Fig. 2 we superimposed an isochrone with 
[Fe/H] =-0.20 (Z=0.012) and an age of 650 Myr. The TO 
fit is very good both in V vs B-V, V vs V-I, and V vs V- 
K diagrams, and the shape of the MS is well reproduced. 
Also the location of #75374 and #76058 stars in the RGB 
is well matched. By this fitting we obtained (m-M)v=14.30 
and E(B-V)=0.43, E(V-I)=0.52, E(V-K)=1.15. 
The reddening value we obtained is compatible with Schlegel 
et al. (1998) maps which, in the direction of NGC 1883, pre- 
dict a E(B-V) value of 0.55, all the way to infinity. 

NGC 1883 turns out to be located 3.9 kpc from the 
Sun towards the anticenter direction. This yields from the 
heliocentric rectangular Galactic coordinates X=-3.71 kpc, 
Y=1.13 kpc, Z=0.41 kpc, and a distance from the Galactic 
center of 12.3 kpc. 

The age we found in this study is somewhat smaller (about 
350 Myr) than in C03, but NGC 1883 can still be considered 
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Table 2. Mean radial velocity, iiiciuborsliij). and parainctois for tlio ol)Scr\'cd stars 



ID 


<RVh> (km/s) 


Membership 


Teff (K) 


log(g) 


vt (km/s) 


[Fe/H] 


79024 


33.9 


NM 










75374 


-31.4 


M 


4172 


2.75 


1.3 


-0.25±0.26 


76058 


-30.3 


M 


4462 


2.75 


1.3 


-0.07±0.46 


81797 


2.8 


NM 










85248 


-18.3 


NM 
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Figure 2. Isochroncs fitting of the photometric data using the 
metallicity obtained in Section 4. left: V vs B-V diagram from 
C03. center: V vs V-I diagram from C03. right: V vs V-K diagram 
(K magnitude from 2MASS). See text for more details. 

one of the few intermediate-age open cluster located in the 
outer Galactic disk. 

The distance estimate is also confirmed by the V mag- 
nitude of the red-clump stars (Girardi & Salaris 2001). 
They are identified with the group of three stars having 
<V>=14.0 and <(B-V)>=1.45 (see Fig. 2). A V magnitude 
of 14.0 for red-clump stars gives a distance of 3.6±0.4 kpc 
from the Sun, that agrees well with the value we obtained 
with isochronc fitting. 

Wc interpret the two member stars at V 13.3 as bright 
RGB staxs. In fact, the interpretation of these two stars as 
clump stars would imply a magnitude difference between 
clump and TO AV^ = 2.1, and, in turn, an age larger than 3 
Gyr (Carraro & Chiosi 19944) . An isochrone for this combi- 
nation of age and metallicity would have too a blue clump, 
and produce a serious mismatch of the evolved stars distri- 
bution. Besides, it would imply a reddening close to zero 
(E(B-V) ^ 0.06), too small for a cluster which in this case 
would lie at 1.7 kpc from the Sun. 



7 SYNTHETIC CMDS 

To better establish the basic parameters of the cluster we 
calculated synthetic V vs B — V CMDs. We started from 
the Girardi ct al. (2000) theoretical isochronc of metallicity 
Z = 0.012 used for our previous fitting. The isochrone was 
shifted in apparent magnitude and colour of the amount 
found in the last Section. The sequence of steps required to 



simulate a CMD aimed to reproduce the NGC 1883 data 
are: 

• The 650 Myr old isochrone oi Z = 0.012 is used to sim- 
ulate a cluster with 3 red clump stars (see previous Section). 
Assuming a Kroupa (2001) IMF, in order to reach this num- 
ber we need an initial cluster mass of about 8.0 x 10^ M©, 
which is assumed herein-after. According to the same IMF, 
we would expect 9-11 stars located beyond the subgiant 
branch region all the way to the asymptotic giant branch 
(AGB) region. The exact number and distribution of these 
stars depend on the simulation seed. Anyhow, this range of 
rmrribor matches the observations. In fact (see Fig. 2), do- 
spite the fact that always some contamination is expected, 
we count 7±3 evolved stars, taking into account Poisson un- 
certainties. 

We have simulated detached binaries, assunung tfiat 30 
percent of the observed objects are binaries with a mass 
ratio between 0.7 and 1.0. This prescription is in agreement 
with several estimates for galactic open clusters (Phelps & 
Janes 1994) and with the observational data for NGC 1818 
and NGC 1866 in the LMC (Elson et al. 1998; Barmina et 
al. 2002). The result of such simulation is shown in Fig. 3(a). 
In this panel we can see that most of the stars - the single 
ones - distribute along the very thin sequence defined by the 
theoretical isochrone. Binaries appear as both (i) a sequence 
of objects roughly parallel to the main sequence of single 
stars, and (ii) some more scattered objects in the evolved 
part of the CMD. 

• In order to estimate the location of foreground and 
background stars, we use the TRILEGAlJ Galaxy model 
code (Girardi et al. 2005). It includes the several Galac- 
tic components - thin and thick disk, halo, and an extinc- 
tion layer - adopting geometric parameters as calibrated by 
Groenewegen et al. (2001) (for other applications see Baume 
et al. 2007, Carraro et al. 2006). The most relevant com- 
ponent in this case is the thin disk, which is modelled by 
exponential density distributions in both vertical and radial 
directions. The radial scale hcigth is kept fixed (2.8 kpc), 
whereas the vertical scale heigth hz increases with the stel- 
lar age t as 

hz =zo{l + t/toT 

with zo = 95 pc, to = 4.4 Gyr, a = 1.66. The simulated 
field has the same area (8.1 x 8.1arcmin^) and galactic co- 
ordinates [i = 163° .08, h = -|-6°.16) as the observed one for 
NGC 1883. The results are shown in Fig. 3(b). It is notewor- 
thy that, in this direction, most of the Galactic field stars 
appear in a sort of diagonal sequence in the CMD. 

t http: / /trilegal.ster. kuleuven.be/cgi-bin/trilegal 
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Figure 3. Synthetic CMD diagram for NGC 1883. Panel a) and 
b): simulated CMDs for the cluster and the field. Panel c) and d): 
simulated CMDs for the cluster and the field with photometric 
errors Panel e): final simulated CMD (cluster+field) . Panel d): 
the observed CMD. 



• Within TRILEGAL, there are various options to take 
reddening into account (Girardi et al. 2005). In the particu- 
lar case of the simulation we are presenting here, we adopt 
an exponential dependence of reddening with distance, cal- 
ibrated at infinity with Schlegel et al. (1998) maps, which 
provide the reddening along a line of sight all the way to 
infinity. 

• We then simulate the photometric errors as a function 

of V magnitude, with typical values derived from our ob- 
servations. The results are shown separately for cluster and 
field stars in panels (c) and (d) of Fig. 3. 

• The sum of field and cluster simulations is shown in 
Fig. 3(e). This can be compared directly to the observed 
data shown in Fig. 3(f). 

The comparison of these two latter panels indicates that 
the selected cluster parameters - age, metallicity, mass, dis- 
tance, reddening, and binary fraction - really lead to an ex- 
cellent description of the observed CMD, when coupled with 
the simulated Galactic field. The most noteworthy aspects 
in this comparison are the location and shape of the turn-off 
and subgiant branch, that are the features most sensitive to 
the cluster age. 

Of course, there are minor discrepancies between the ob- 
served and simulated data, namely: (i) The simulated cluster 
is better delineated in the CMD than the data. This may 
be ascribed to a possible underestimate of the photometric 
errors in our simulations, and to the possible presence of dif- 
ferential reddening across the cluster (ii) There is a deficit 



of simulated field stars, that can be noticed more clearly for 
V < 16 and {B — V) < 1. This is caused by the simpli- 
fied way in which the thin disk is included in the Galactic 
model: it is represented by means of simple exponentially- 
decreasing stellar densities in both radial and vertical direc- 
tions, and does not include features such as spiral arms, 
intervening clusters, etc., that axe necessary to correctly 
describe fields at low galactic latitudes. Anyway, the fore- 
ground/background simulation we present is only meant to 
give us an idea of the expected location of field stars in the 
CMD. 

Although these shortcomings in our simulations might 

probably be eliminated with the use of slightly different 
prescriptions, they do not affect our main results, that re- 
gard the choice of cluster parameters Then, we conclude 
that (m-M)v = 14.3, Eb-v = 0.43, 650 Myr, Z = 0.012 
([Fe/H] = —0.20) well represent the cluster parameters. All 
these values are uncertain to some extent: 

• From isochrones fitting we can estimate a maximum 
error of 10 percent (70 Myr) in the age. One should keep 

in mind, however, that the absolute age value we derived 
depends on the choice of evolutionary models, and specially 
on the prescription for the extent of convective cores. For 
the stellar masses involved (A/to ~ 1.8 Mq for NGC 1883), 
our models (Girardi et al. 2000) include a moderate amount 
of core overshooting. 

• Our best fit model corresponds to E(B-V)=0.43 and 
(m-M)\/=14.30, but typical errors for both measurements 
can reach 10-15%. 

• The initial mass estimate depends heavily on the choice 
of IMF, that determines the mass fraction locked into low- 
mass (unobserved) objects. The value of 8.0 x 10^ Mq was 
obtained with a Kroupa (2001) IMF, corrected in the lowest 
mass interval according to Chabricr (2001; details arc given 
in Groenewegen et al. 2001), and should be considered just 
as a first guess. At present ages, supernovae explosions, stel- 
lar mass loss, and tidal efi^ects (Lamers et al. 2005) would 
have reduced this mass by about 40 percent. 



8 SURFACE DENSITY PROFILE 

According to C03 and Dias et al. (2002), NGC 1883 has a 
radius of 2.5-3.0 axcmin. This is the distance from the clus- 
ter center where star counts remain flat, reaching the level of 
the Galactic field. However, to provide more physical quanti- 
ties, we perform here a comparison with King (1962) models 
(Fig. 4). 

To this aim we first adopt a new cluster center based on 
the stellar concentration. C03 chose star #12 as center (see 
Fig. 5), but a more detailed inspection of the cluster appear- 
ance suggests that a more appropriate option for the center 
is X = 495, Y = 605 pixels (see Fig. 5), which in turns is 
(a20oo = 05''25'"58''; Jaooo = +46''29'03"). 
By adopting the new cluster center, we counted the num- 
ber of stars as a function of the distance from the center, 
and correct these counts by background contamination. The 
choice of the background is illustrated in Fig. 5, and implies 
that one can expect to encounter ~ 4.5 stars/arcmin^ in the 
direction of NGC 1883 down to V^20. This value is in rea- 
sonable agreement with the predictions of the TRILEGAL 
model (Section 7), according to which we expect to find 5.2 
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Figure 4. Surface density profile for the cluster. The solid line is 
a King model for the indicated parameters. T,g shows the level 
of the field background. 



stars/arcmin^ towards at (l=163.1,b=+6.2). 
In Fig. 4 we show star counts together with their Poisson un- 
certainties. The sohd hne is the best-fit King profile, which 
yields a core radius rc ~ 0.30 arcmin and a concentration 
c — log{^) = 1.43, where rt is the cluster tidal radius. 
rt results 8.1 arcmin, much larger than the region covered 
by our photometry. These values can be compared with the 
recent study by Piskunov et al. (2007), where the authors 
study the rc and rt distribution of a sample of open clusters 
in the Sun vicinity. According to this study, the open clus- 
ters mean c value is ~0.6. Therefore, NGC 1883 looks like a 
poorly concentrated cluster, on the verge of dissolving into 
the Galactic disk general field. 
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Figure 5. The field used for the Cluster Luminosity Function 
determination. The Cluster Region and the Comparison Field 
are indicated. 



Table 3. Completeness results in NGC 1883. 



V 



% 



10.5 
11.5 
12.5 
13.5 
14.5 
15.5 
16.5 
17.5 
18.5 
19.5 
20.5 



100.0 
100.0 
99.9 
99.5 
97.3 
96.1 
95.8 
90.1 
73.4 
62.0 
41.3 



9 CLUSTER LUMINOSITY AND MASS 
FUNCTION 

With the aim of deriving an estimate of NGC 1883 inte- 
grated luminosity and mass, we construct the luminosity 
function (LF) of the cluster. 

As a first step, we considered the Cluster Region and the 
Comparson Field: the first one is a circle centered in the 
adopted cluster center and 2.'5 in radius; and the later is a 
region limited for a circle of 3.'ll in radius and a 7.'l x 7.'l 
square (see Fig. 5). Both regions were chosen to cover the 
same area. 

We estimated then the completeness of our data as in 
Baume et al. (2004) and assume that it is the same across 
the whole region. This seems to be a good approximation, 
provided the cluster is not very crowded. The completeness 
as a function of magnitude V is given in Table 3 

Besides, we assume that our comparison field provides a 
good estimate of the contamination by field interlopers and 
is valid across all the cluster surface. The results are plotted 



in Fig. 6. 

Here, the completeness, background corrected stars counts 
as a function of magnitude are shown together with their 
Poisson uncertainties. The LF is a rising function over the 
magnitude range for which our completeness is larger than 
50%, except for the last bin, at V larger than 18. This seems 
to suggest that evaporation of low mass stars has been im- 
portant during NGC 1883 evolution, in agreement with the 
findings in Sect. 8. Since this happens at the limit of our 
photometry, a deeper study is mandatory to make stronger 
our results. 

We note in addition the presence of a statistically significant 
depression at V 15.5, where also the MS in the CMD (see 
Fig. 1) does show a gap. Such gaps are common in open 
clusters, and a variety of explanations have been discussed 
over the years for their occurrence (Rachford & Canterna 
2000). 

Apparent magnitudes were then transformed into the abso- 
lute ones using the cluster apparent distance modulus (m- 
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Figure 6. Cluster Luminosity Function. 
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Figure 7. Cluster Mass Functions. 



M)y = 14.3. This allows us to compute the integrated abso- 
lute magnitude of the cluster by adding up star fluxes. We 
obtain Mv^ = -2.4±0.2. This value has to be considered as 
an upper limit, due to the possible presence of several field 
stars which we are not able to completely remove. 



Figure 8. Galactic abundance gradient. Open squares are clus- 
ters from Friel et al. (2002), while open triangles are clusters 
found in the WEBDA database (see text for more details). The 
open star is the value for NGC 1883. The position of the Sun 
(8.5,0) is indicated 



present day MF is shown in Fig. 7. A linear fit through the 
points yields 2.1±0.5 for the a index, close to the Salpeter 
2.35 value. As for the LF, we notice a deficiency of less mas- 
sive stars, presumably due tidal mass loss. The integration of 
the present day MF gives a current total mass for NGC 1883 
of about 240 M© , 3 times smaller than the mass at birth we 
derived in Sect. 7. This lend further support to a scenario in 
which substantial mass loss occurred as the cluster orbited 
the Milky Way. 



10 DISCUSSION AND CONCLUSIONS 

In Figure 8 we plot the Galactic radial abundance gradient 
as derived from Friel et al. (2002, hereafter F02) from spec- 
troscopic observations of a sample of open clusters (open 

squares). Wc added also a sample of clusters not observed 
by F02 (open triangles). We selected these objects from the 

WEBDA database^ with the condition of having metallici- 
ties spectroscopically dotorminatod to be homogeneous with 
F02 and with this paper. We plotted only clusters having 
age comparable with NGC 1883 (0.5<Age<2.0 Gyr, the 
same range used in F02 to identify the intermediate-age clus- 
ters sample). The whole sample defines an overall slope of 
-0.06±0.01 dcx/kpc (dashed line). 

The open star represents NGC 1883, which follows the 
general trend very nicely. In fact, the continuous line, which 



Finally the My distribution was converted into mass func- 
tion (MF) using Girardi et al. (2001) theoretical models. The 



8 WEBDA database on open clusters can be found at 
http://www.univie.ac.at/webda/ 
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represents the radial abundance gradient determined by in- 
cluding also NGC 1883 basically coincides with the previous 
one, having a slope of -0.05±0.01 dex/kpc. But we must no- 
tice that the slope is mainly driven by the two inner most 
metal rich open clusters (the inner triangles) . If wc exclude 
these points, the gradient is flat having a value of -O.OliO.Ol 
dex/kpc (the dotted line). Considering also the dispersion 
of the data and the error of the fit, we are led to the conclu- 
sion that Galactic gradient is flat, at least outside the solar 
circle. Inside it probably grows up to supor-solar mctallic- 
ity ([Fe/H]'~-|-0.2). However this last point is not sure being 
based only on the metallicities of two clusters, and it needs 
further investigations. 

On the other hand, as discussed by F02, the oldest open 
clusters (Age>4.0 Gr) formed in our Galaxy demonstrate a 
strong gradient also outside the solar circle, as shown in Fig. 
3 of F02. This is compatible with the scenario of a primor- 
dial pollution of the interstellar matter from which clusters 
were formed mainly in the inner regions of the Galaxy, fol- 
lowing which the chemical enrichment was extended to the 
outer regions. The results was a present day homogeniza- 
tion of the chemical content in the Galactic disk giving a 
flat gradient as defined by the younger open clusters outside 
the solar circle. However this scenario is not valid for the 
outer regions of the Milky Way. 

Carraro et al. (2004, hereafter C04), and Yong et al. (2005) 
found the presence in the outer disk of old open clusters 
(Age>4 Gyr) having a mctallicity of about -0.5 dex, well 
above the value of the gradient as defined by the inner old 
clusters. This means that probably the outer disk had a pe- 
culiar chemical evolution, probably driven by the contami- 
nation of external objects (dwarf galeixies) cannibalized by 
the Milky. 

On the other hand only few intermediate-age clusters 
were observed for metallicity determination in the outer disk 
at medium-high resolution, the majority of them being lo- 
cated close to the Sun. 

The addition of NGC 1883 indicates that the gradients 
keep flat for half a kpc more, suggesting a trend which fur- 
ther studies can possibly conflrm. 

We finally note that the gradient exhibits quite a signif- 
icant scatter. One may wonder whether this solely depends 
on observational errors, or whether this scatter reflects a 
true chemical inhomogeneity in the Galactic disk. 

Further observations are needed, especially in the not 
well explored Galactic region between 14 and 22 kpc from 
the Galactic center. We are actually engaged in such explo- 
ration, with new results being presented in Carraro et al. 
(2007, submitted to AJ). 
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